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SUMMARY 

Performance of a 30 GHz FDMA uplink to a processing satellite is 
modelled for the case where the on-board demultiplexer is implemented 
optically. Included in the perforance model are the effects of adja- 
cent channel interference, intersymbol interference, and spurious signals 
associated with the optical implementation. Demultiplexer parameters are 
optimized to provide the minimum bit error probability at a given 
bandwidth efficiency when filtered QPSK modulation is employed. 


INTRODUCTION 

Satellite communication using frequency division multiple access 
(FDMA) on the uplinks and time division multiple access (TDMA) on the 
downlinks has attracted much interest (1,2). FDMA on the uplink permits 
the use of ground transmitters that do not require amplifiers having 
excessively high power. Also, FDMA does not require complicated network 
timing. TDMA on the downlink takes advantage of recent developments in 
satellite on-board processing and switching capabilities to provide high 
data rate downlinks to VSAT-type ground receivers. In addition, the 
heavily-used C-band and Ku-band frequencies will be supplemented by 
higher frequency Ka-band transmission (30 GHz uplink / 20 GHz downlink). 
This permits the use of smaller ground terminal antennas, but at a cost 
of higher rain attenuation. 

On-board processing is needed to efficiently service multiple users 
while at the same time minimizing earth station complexity. Figure 1 is 
a simplified overview of a SATCOM system that services FDMA uplink users. 
The processing satellite first receives the wideband uplink at 30 GHz and 
downcon-verts it to a suitable IF. A demultiplexer then separates the 
composite IF signal into assigned channels. All channels are then de- 
modulated by "bulk" demodulators, with the baseband signals being routed 
to the downlink processor for retransmission to the receiving earth sta 
tions via a high-rate TDMA 20 GHz downlink. This type of processing 
circumvents many of the difficulties associated with bent-pipe repeaters. 
First, uplink signal distortion and interference are not retransmitted on 
the downlink. Second, downlink power can be allocated in accordance with 
user needs, independent of uplink transmissions. This allows the uplink 
users to employ different data rates as well as different modulation and 
coding schemes. In addition, all downlink users will then have a common 
frequency standard and symbol clock on the satellite, which is useful for 
•network synchronization. 
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These considerations led to a requirement for on-board multi- 
channel demodulators (MCD) that can separate and process the individual 
transmissions with minimal degradation in bit error probability. 
Implementation of an MCD is critical because future systems will be 
highly bandwidth-efficient, which implies very close spacing of the 
carriers in the composite FDMA uplink. 

On-board FDMA demultiplexers can be implemented in a variety of 
ways. One way is to do a wideband A/D conversion on the uplink signal 
received at the satellite, followed by digital processing that performs 
the channel filtering and demodulation operations [3]. However, on-board 
demultiplexing can also be performed using integrated optics [4,51. An 
acousto-optical spectrum analyzer performs both down-conversion" and 
c annel filtering, with potential savings in hardware size and weight. 

This paper shows how an acousto-optical demultiplexer can be 
modelled in system performance analyses. Bit error performance is deter- 
mined in the presence of adjacent channel interference, intersymbol 
interference, and spurious signals generated by the optical processing. 

ON-BOARD DEMULTIPLEXER 

. . . An acousto- °Ptical spectrum analyzer (Fig. 2) employing heterodyne 
detection can function as a channelized receiver. The spectrum analyzer 
converts the composite FDMA uplink into acoustic waves in a Bragg cell. 
These acoustic waves modulate a laser beam, and diffract the beam at 
ang es proportional to the uplink RF signal frequencies. Reference beams 
are also provided to achieve heterodyne operation, resulting in larqer 
dynamic range. The diffracted light impinges on an array of photodetec- 
tors, which function as square-law detectors, and the individual photo- 
currents are routed to QPSK demodulators. Thus, the acousto-optical 
spectrum analyzer serves as both a channelizer and downconverter , so that 
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Fiq. 2 Acousto-opt led 1 Heterodyne Spectrum Analyzer 


the composite uplink signal is demultiplexed into separate channels, each 
at a common IF. 

To estimate system performance from a demultiplexer of this type, 
it is first necessary to determine its transfer function. In a classical 
linear system, a sinusoidal input to the system results in a sinusoidal 
output at the same frequency, whose amplitude and phase depend on the 
frequency. But for the heterodyne system considered here, a sinusoidal 
input results in not only a sinusoidal output at that frequency, but also 
sinusoids at f + ^ nF, where F is the channel spacing. These spurious 
sinusoids are generated internal to the demultiplexer by the reference 
frequency comb. In contrast, the frequencies provided by other transmit- 
ters external to the demultiplexer make up the ACI, which is characteris- 
tic of all FDMA systems. 

DEMULTIPLEXER TRANSFER FUNCTION 

As indicated in Fig. 2, two channels constitute the acousto-opt ical 
spectrum analyzer: the "signal" channel and the "reference" channel. 
When the beam has a Gaussian cross section, the light into the signal 
channel Bragg cell can be expressed as 


e -csP+>Wi. i (1) 
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where ? L is the light frequency, x denotes distance from the center of 
the Bragg cell, and c ff is a constant determined by the laser beamwidth. 
This light is modulated by an acoustic wave produced by the input 
sinusoid of frequency f s , which can be expressed as * 

e'WsO-J) (2) 

where v is the acoustic velocity in the Bragg cell. Therefore, the modu- 
lated light out of the cell is the product of (1) and (2), and the signal 

channel light distribution in the k-plane is the spatial Fourier trans- 
form: 

r s (k) = e , 2 "<h+/s)< r (3) 


where d s is the length of the Bragg cell. Similarly, the light distribu- 
tion resulting from the multi-diffracted beam in the reference channel is 

00 ,±a 

r R (k )= Y, c' 29 UL +t*"' F " l € -c R *’- t 2* { i*±2L+ k] Tdx * 

n = - oo * — 

One of the reference beams is directed toward the signal channel light 
distribution (3) in the k-plane. However, because the other reference 
beams also overlap the signal beam to some extent, spurious output 
signals occur. If a high bandwidth efficiency is required, then the 
beams must overlap more, which implies a higher level of spurious 
signals. 

The total light intensity in the k-plane is 


I Fs + F r \ 2 = \F S \ 2 + \F r \ 2 + 2 ReF s F' R (5) 

By suitably filtering the photodetector output, the only important 
contribution to the output will be the cross product G(k)= Re F F* . The 
photocurrent is proportional to the integral of the intensity over the 
photosensitive area: 


I = 


L 


*0--? 


*o-f 


G(k)dk 


( 6 ) 


S Sf tL th p e hoto C =en° n °Nof e 1 e P t h 0 f tOCe b 1 e 1 th" * nd K is the 

frequency and let f = f f iet f ® nominal channel carrier 

this nominal frequency. \ssume That The k^laT^i fr< L < * uenc > r relative to 
to the frequency f A Y so that k = , *" pl *£ e locatlon k o corresponds 

translated in frequency to baseband, reduces to" P hot ^®ll output. 


oo 

7 = Hft(f> n F)cos2x(f-nF)t 

n= — oo 


( 7 ) 


182 



where 


>1 / ) 


Hji 

'• I • r 

[ ? n m \\ n i\ { r 

HIS — (j X 4- ri/'v) — 


/< 1 

f ns - ( / r 

l 

MM n /\'( z -j y ) 

n! >j) 

I 4 ij 

1 V 

<ir d y 




rn 


The function H^(f,0) is the transfer function winch will be denoted f ) , 
and which can be calculated by numerical integration after the system pa- 
rameters have been selected. The terms for which n = 0 give the amplitudes 
of the spurious sinusoids. 

Figure 4 is the computed transfer function for a particular set of 
MCD parameters. The parameters were selected to give the minimum bit 
error probability in the presence of ACI and ISI when the earth station 
transmissions are filtered with fifth-order Butterworth filters having 
time-bandwidth products of 0.5. A bandwidth efficiency of 1.6 bps/Hz has 
been assumed. We now describe the method used to compute the ACI and 
ISI, and then the bit error probability itself. 




riq. 1 Interference Effects 
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ADJACENT CHANNEL INTERFERENCE 

Design of any bandwidth-efficient FDMA system involves a fundamen- 
tal trade-off. If the system bandwidth is narrow, we achieve good ACI 
performance at the cost of high ISI. Widening the bandwidth reduces the 
ISI but increases the ACI. The design procedure is generally to select 
filter types and bandwidths that give the best bit error performance in 
the presence of both ACI and ISI. We first consider the ACI. 

Each ground transmitter is assumed to include bandpass. filtering to 
reduce the amount of ACI entering the satellite receiver. As a first 
approximation for performance analysis, the demultiplexer can be treated 
conventionally as a bank of bandpass filters, each followed by a demodu- 
lator. The model is later generalized to take into account the spurious 
signals that are characteristic of the demultiplexer implementation. 

It is straightforward to compute the ACI under the assumption that 
the interference can be treated as noise that adds to the thermal noise 
at the receiver input. This assumption is valid when there is a large 
number of interfering users because according to the central limit 
theorem, this implies that the interference has nearly Gaussian statis- 
tics. We also assume for simplicity that all transmissions arrive at the 
satellite with equal power. 
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S(f - »F) 


sin 2i rT(/ - nF) 
2 *T(f - nfj . 


( 10 ) 


All around transmitters have identical filters, and the n th filter 
Ail grounu Thus the filtered transmission 

transfer function is denoted by H T (f-nF). Thus cne iniu _ f f-nFh z 

from the n-th interferer has a spectral density given by S(fnF)|l T U- H) 
Sunnose the transfer function of the on-board demultiplexer is H^(f), 
which was evaluated in the previous section. Then the spectral density 
of the interference into the demodulator is 


X>'(/ - nF)\H T (f - nF)Il R (f)\ 2 


Assume that the symbol detector is a f liter matched to th 
undistorted symbol (i.e. an integrate-and-dump _ detector) .Its tran 
function is H , ( f ) =sin (2*Tf ) / 2t!Tf. Then the total ACI power out ol the 
matched f i lte? F relative to the undistorted signal power is 


/ oo 

S(f - nF)\H T (f - nF)n R (f)H MF (f)\ 2 df 
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This is added to the thermal noise to estimate the error probability. 


INTERSYMBOL INTERFERENCE 

Unlike the ACI, the ISI cannot be accurately approximated as addi- 
tive Gaussian noise. Instead, we determine explicitly the effect of the 
transmitter and receiver filtering on the amplitude of the signal out of 
the integrate-and-dump detector. 

Expressing any one of the unfiltered QPSK signals before transmis- 
sion as s ( t ) = m (t) cos ( 2 f t + ), the complex modulation is 
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where a and b are binary data on the I- and Q- channels respectively, 
which are assumed to be offset l/2-symbol. The modulation spectrum of 
the filtered and demultiplexed signal at the input to its demodulator is 
fi??) = M (f)H Jf)H.(£), where M(f) is the spectrum of the undistorted QPSK 
modulation. The Fourier transform of this is the distorted modulati, 
which we will call m(t), and the output of the I-channel integrate and 
dump detector in the demodulator is 

L r *(0* - r £>- £ *W«/> (l^) 1 <“> 
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alter changing the integration variable. This shows th^t- +-k q ^ *. 
output includes contributions from not only the desired (n-0) symbol* 0 ? 0 ) 
from all symbols which is what is mean/ by intersymboi /nte'rf erence 
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By differentiation, we find that the value of for which tho 
average value of S is a maximum is the solution of 
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Once the optimum 'f is found we can compute the contribution to S from tne 
n-0 symbol and from all important interfering symbols. It has been found 
that the n - 1,2, 3, and -1 interfering symbols are the important ones in 
our application. 

The bit error probability is then computed by averaging the error 
probability , conditioned on a particular sequence of these interfering 
symbols, over all 16 possible sequences of these symbols. The noise in 
this computation consists of the ACI, which was calculated in the pre- 
vious section, plus the thermal noise at the input to the satellite 
receiver. If N 0 is the thermal noise density at the input, the thermal 
noise power at the output of the integrate— and— dump detector is 


r 00 

^0 / \H R (!)H MF {f)\ 1 df 

2 — oo 


(17) 


which is added to the ACI in the error probability computation. 


RESULTS 


Having derived the 
calculated the ACI and 
expression for the total 


demultiplexer transfer function, and having 
ISI, it is .straightforward to write down an 
power out of the device: 


P =HT, [°° S U ~nF)\R T {J -nF)H R U,mF)n M F{f)\ 7 d[ (18) 
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The first sum is over the number of interfering earth stations, while the 
second is over the number of spurious signals produced by the reference 
beams. The m=n=0 term is the desired output; the m-0, nfO terms are the 
ACI* the m^O n=0 terms are the spurious outputs when there is only one 
uplink signal; and the m* 0 , njo terms are additional spurious - outputs 
arising from interaction of the interfering channels with the n^O refer- 
ence beams. Of this latter category, the m=n terms dominate, so they are 
included in the performance calculation. The mfn interaction terms are 

illustrated^in Frg.^2- t he ^ GrrQr pro bability when the MCD transfer 

function is that shown in Fig. 4- Fifth-order Butterworth filters are 
used in the earth station transmitters, and integrate and dump symbol 
detection is assumed. The effects of ACI, spurious signals, and ISI are 
all included in the calculation. Also shown for comparison is a case 
where the optical MCD is replaced by a bank of bandpass fifth order 
Butterworth filters whose bandwidths are optimized to give the smallest 
bit error probability in the presence of the same interference Figure 
6 shows that the performance of an optical MCDi compares favora y wi 
that of an MCD implemented electronically. 

Figure 7 shows the sensitivity of system performance to timing 
errors in the symbol detection circuitry. This is evaluated by varying 
the sampling time about its optimum value when computing the error pro 
bability. It is evident that when the timing error is less than approxi 
mately ten percent of the symbol duration, its contribution to bit error 
degradation is insignificant compared to the ACI, spurious signals, and 

ISI. 



Fiq. 5 Analytical Model for Performance Evaluation 
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Fig. 6 Bit Error Probability 


Eig. 7 Sensitivity to Symbol Synchronization Error 


DISCUSSION AND CONCLUSION 

The goal of this effort was to evaluate the performance of an 
optically implemented on-board demultiplexer that can service inex- 
pensive, low power earth stations. For such earth stations, no attempt 

V im P rove bandwidth efficiency by using advanced modu- 
lation and coding schemes. in addition, we have assumed that no effort 
has been made to minimize ISI by careful design of the earth station 
transmitter filters. We have shown that an optically implemented MCD 
which promises size and weight advantages over other implementations of 
on board processors, can perform as well as other implementations so far 
as bit error probability is concerned. 

However, to- achieve comparable performance with the optical MCD 
whose transfer function must be carefully controlled, it was found 
necessary to reduce the ground transmitter filter bandwidths. As seen in 
Fig. 6 , this worsens performance at low signal-to-noise ratios because of 

^ e . rea ff5 P 1 t 5 at results - The optical MCD transfer function (Fig. 
4 ) has amplitude and phase characteristics different from most classical 

r mK P ° n J eS ’ *.u 1S °' s P urious responses from the reference beams 
There J ore ' t he ACI/ISI trade-off for an optical implementation 
differs significantly from that for an electronic one. 
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